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Abstract 

The fluorescence yield data used for shower reconstruction in the Auger, 
HiRes and TA experiments are different, not only in the overall absolute 
value but also in the wavelength spectrum and the various atmospheric de- 
pendencies. The effect on the energy reconstruction of using different flu- 
orescence yield parameterizations is discussed. In addition, the impact of 
a change in the fluorescence spectrum depends on the optical efficiency 
of the telescopes. A simple analytical procedure allows us to evaluate the 
combined effect of fluorescence yield and optical efficiency showing a non- 
negligible deviation between the energy scales of TA and Auger. However 
no relevant effect is found in the comparison between HiRes and Auger. Fi- 
nally we show that a similar procedure could also be applied with real data. 



1 Introduction 



The accurate measurement of the calorimetric energy of UHECR showers by 
means of fluorescence telescopes relies on the precise knowledge of the air-fluorescence 
yield. In practice, the fluorescence yield is described by a set of parameters that al- 
lows us to convert the fluorescence intensity recorded by the telescope into energy 
deposited in the atmosphere at given atmospheric conditions. Auger, HiRes and 
Telescope Array (TA) use different parameterizations of the fluorescence yield 
(absolute value, wavelength spectrum and atmospheric dependency). The com- 
parison of the energy spectra of Auger and HiRes/TA indicates a possible dis- 
agreement in the energy scale 0], which could be partly due to using different 
fluorescence yield data. In this work we evaluate the effect of the fluorescence 
yield choice on the relative energy scale of these experiments. 

While deviations in the absolute value of the fluorescence yield translate di- 
rectly into deviations in the reconstructed energy, the effect of differences in the 
relative wavelength spectra depends on the optical efficiency of the telescope. 
Therefore, in order to evaluate the actual influence of the fluorescence yield on 
the energy scales of these experiments, the combined effect of the assumed Y A 
values (i.e., absolute spectral values) and the optical efficiency s A has to be anal- 
ysed. 

The dependence of the fluorescence yield with the atmospheric parameters is 
determined by the characteristic pressures P' r Neglecting the humidity contribu- 
tion to quenching and the temperature dependence of the collisional cross section, 
the differences in the assumed P' A values do not lead to significant disagreements 
in the energy reconstruction, since most of the shower energy is deposited in the 
lower atmospheric layers where quenching effects saturate. However these hu- 
midity and temperature contributions to P' A give rise to a non-negligible effect on 
the reconstructed energy [0. 

The wavelength spectrum assumed by Auger, HiRes and TA, as well as the 
corresponding optical efficiencies of these experiments are described in section [2j 
The analytical method we used to compute the differences in the energy scales 
is described in section [3] and the corresponding results are presented in section |4j 
This theoretical method overestimate somewhat the energy deviations with respect 
to those from a full reconstruction of real showers. In section [5] we discuss on the 
applicability of this procedure to real data. 
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Figure 1: Absolute fluorescence yield (dry air) spectra at 1013 hPa and 293 K employed by the 
Auger, HiRes and TA experiments (blue bars). The red line represents the optical efficiency of the 
corresponding telescopes including all optical elements 

2 Fluorescence Yield datasets and optical efficien- 
cies 

Data on the fluorescence yield used by Auger, HiRes and TA and the optical ef- 
ficiency of the corresponding telescopes are decried below. These parameters are 
shown together in Figure[T]for comparison. 

2.1 Auger 

The fluorescence yield model presently used by the Pierre Auger Collaboration is 
defined as follows. The absolute value of the yield for the 337 nm band Y 331 is 
the one measured by Nagano et al. ||3). The yields for the remaining wavelengths 
are distributed according to the relative intensities measured by the AIRFLY Col- 
laboration 01. The characteristic pressures P' A are also those from AIRFLY flU. 
Auger is the only experiment taking into account the water vapor quenching of the 
atmosphere and the temperature dependence of the collisional cross section ||5). 

The relative optical efficiency of the Auger telescopes O is dominated by the 
filter (MUG6 glass manufactured by Schott-Desag) although other elements (i.e., 
mirror, corrector ring, and PMT) also contribute. 

2.2 HiRes - TA (Middle Drum) 

The HiRes Collaboration uses the absolute value of the integrated fluorescence 
yield (300 - 400 nm) from Kakimoto et al. Q. The relative intensities of the 
337, 351 and 391 nm bands are those measured by Kakimoto while the intensi- 
ties of the remaining spectral components are distributed according to those from 
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Bunner ||S). The P' A values are those reported by Kakimoto 0. 

The optical efficiency of the HiRes telescopes is described in [9]. The filter is 
very similar to the one used by Auger. One of the fluorescence detectors of HiRes 
has been refurbished and is presently used by the Telescope Array Collaboration, 
i.e., the so-called Middle Drum (MD) telescope. 

2.3 Telescope Array (Black Rock Mesa & Long Ridge) 

The Black Rock Mesa (BRM) and Long Ridge (LR) eyes of the TA Collaboration 
use the absolute value of the integrated fluorescence yield (300 - 400 nm ) from 
Kakimoto et al. [7], that is, the same absolute reference of HiRes. However TA 
uses the wavelength spectrum measured by FLASH [10]. This spectrum has two 
additional bands not reported by Kakimoto that are added by TA to get their total 
absolute yield. The P' A values are also those from Kakimoto et al. 0. 

These BRM and LR telescopes use a BG3 glass (manufactured by Schott) as 
UV filter. This filter has a higher transmittance over a wider wavelength interval 
than the MUG-6 glass employed by Auger. This is a relevant feature that poten- 
tially leads to a lower signal-to-noise ratio. In addition, other molecular bands, 
not taken into account because they are absorbed by other narrower UV filters, 
can contribute to the total fluorescence signal recorded by TA. 

3 Method 

A first estimate of the effect of a change in the fluorescence yield on the re- 
constructed energy can be obtained as follows. Neglecting other contributions 
(e.g., atmospheric transmission and geometrical factors), the fluorescence sig- 
nal recorded by a telescope due to a true energy deposition E, is proportional 
to E, Yia Yai^a, where s A is the optical efficiency of the telescope and is the 
true fluorescence yield, both dependent on wavelength A. Hereafter we relax the 
notation using Ye = Yja Ya£a- 

Let us assume two experiments with efficiencies s x and s 2 using different flu- 
orescence yields Yi and Y 2 respectively. For a given E t value, the first experiment 
reconstructs an energy E\. If the fluorescence yield is replaced by Y 2 , the recon- 
structed energy would be E 2 i, that is, 

E t Y t s l = E l Y 1 s l = E 21 Y 2 s 1 . (1) 
Therefore the effect of a change in the fluorescence yield is given by 
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£21 = Y^ 

Ex Y 2fil V 
A relationship similar to ([T]) applies to experiment 2 

E t Y t £ 2 = £ 2 Y 2 £ 2 . (3) 

In principle the ratio of the energies measured by these two telescopes E 2 /E\ 
cannot be extracted from the above expressions. Nevertheless if the true fluores- 
cence yield is known, expressions ([T]) and ([3]) would lead to the solution 

£2 = Y1S1 Y t e 2 
£1 Y 2 e 2 Yjfi! ' 

and thus, in case one of the two experiments is assumed to use the true fluores- 
cence yield (e.g., Yi = Y t ), equation ([4]) will give the simple result 

^1 = ^1. (5 ) 
Ei Y 2 e 2 

This is similar to equation ([2]), although for the combined effect of efficiency and 
yield the ratio has to be computed with the efficiency of the experiment using the 
wrong yield. Note that even if the assumed true yield turns out to be wrong by a 
scale factor, expressions ([4]) and ([5]) are still valid, that is Y t only needs to describe 
the true relative fluorescence spectrum. 

The above procedure allows the evaluation of the deviation in the energy de- 
position at a given atmospheric depth due to a change in the fluorescence yield 
([2]) or in both fluorescence yield and optical efficiency Q[5]). The corresponding 
effect on the total calorimetric energy of the shower can also be calculated from 
the above expressions as long as we replace the above dentition, Ye = Yja Y a s a , 
by the following one 

dX^G(X) Y a {X)s a T a {X) , (6) 

where dn/dX is the normalised development of energy deposited in the atmo- 
sphere by the shower, G is the geometrical factor, X is the slant depth and T A is 
the atmospheric transmission of the light traversing the path from the particular 
shower point to the telescope location. Using ([6]) has allowed us to calculate more 
realistically the energy deviations. For the development of deposited energy, a 
Gaisser-Hillas profile has been used and the T A values have been computed from 
ifTTTl (see E0 for details). 
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Table 1: Energy deviation with respect to Auger AE, defined in ([7]), for differ- 
ent optical efficiencies (rows) and fluorescence yields (columns). Values (black) 
in the upper row represent the energy deviation if Auger would have used the 
fluorescence model of TA or HiRes. Assuming that Auger uses the true fluores- 
cence spectrum, the energy scale of HiRes (TA) is shifted with respect to that of 
Auger by the blue (red) value. Results are shown for three cases: fixed pressure 
and temperature, realistic Gaisser-Hillas shower without and with temperature and 
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4 Results 



In the first place we present the results neglecting Rayleigh scattering and geomet- 
rical factors. Using the data of section 2, we have obtained the Ye values for all 
the yield-efficiency combinations. The percentage energy deviations with respect 
to that of Auger, defined in this work as 



AE = 200 x 



£ _ £ Auger 



(7) 



E + £Auger ' 

have been evaluated at fixed atmospheric conditions (800 hPa and 293 K) for both 
HiRes and TA (see table 1). The results in the first row have been obtained using 
the optical efficiency of Auger and therefore they represent the energy effect if 
the Auger yield were replaced by those of the other experiments ([2]). Those of 
second and third rows have been evaluated using the efficiencies of HiRes and TA 
respectively and therefore they give us the relative energy scale of these experi- 
ments with respect to Auger if we assume that Auger is using the true fluorescence 
spectrum ([5]). 

Next a more realistic calculation has been carried out by averaging the energy 
ratios on the shower development including the Rayleigh effect and geometrical 
factors, that is, we have used the Ye definition given in expression ([6]). Calcula- 
tions for a vertical 10 19 eV proton shower as a function of the distance between 
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Figure 2: Relative deviations of the energy scales as a function of the distance between the 
telescope and the shower maximum for a vertical proton shower of 10 19 eV. Rayleigh scattering 
has been included in this calculation. 



the telescope and the shower maximum in the range 600 - 2600 g cm" 2 have been 
performed. The results for the above yield-efficiency combinations are plotted in 
figure [2j As shown in the figure, AE depends on the mass thickness between the 
telescope and the shower maximum due to the Rayleigh effect. The corresponding 
AE values for a typical vertical shower at 15km from the telescope are also shown 
in the table (columns 3 and 4). 

These calculations have been repeated including the temperature and humidity 
contributions to the P' A values in Auger with typical atmospheric profiles (see Q 
for details). As shown in the table (columns 5 and 6), the disagreement between 
the relative energy scales for both cases HiRes-Auger and TA-Auger is partly 
cancelled. This is expected since the reconstructed energy increases when the 
temperature and humidity is included. 

In summary, assuming that the AIRFLY fluorescence spectrum is the true one, 
according to our calculations the relative energy scale of TA is 13% (10%) larger 
than that of Auger if the temperature and humidity effects are (not) included in 
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the Auger yield. On the other hand, the Auger energy would increase by 10% 
(6%) if this experiment replaced its yield by that of TA, without (with) including 
temperature and humidity effects. Finally we infer that the energy scale of TA is 
shifted by about 10% with respect to that of HiRes. 

5 Towards a precise determination of the relative 
energy scales 

The above results on the relative energy scales of Auger, HiRes and TA have been 
obtained assuming ideal showers that only emit UV radiation from atmospheric 
fluorescence. In practice a non-negligible amount of Cherenkov light is emitted 
and partly detected by the telescope. The contribution of the Cherenkov light has 
to be included in the reconstruction algorithm to measure the shower energy [fT2|. 
The larger the Cherenkov contribution, the smaller the effect of a change in the 
fluorescence yield on the reconstructed calorimetric energy. As a consequence, 
a given change in the fluorescence yield does not translate entirely in the energy 
scale and therefore the above results are expected to overestimate the energy effect 
in real data. 

The effect of a change in the fluorescence yield can easily be obtained from the 
observed deviation in the reconstructed energy of real showers when replacing in 
the reconstruction algorithm the Y A data by that of other experiment. However, the 
evaluation of the combined effect of optical efficiency and fluorescence yield is not 
trivial since replacing the optical efficiency in the reconstruction algorithm by that 
of a different experiment gives rise to a wrong reconstructed energy. However, 
as shown next, the problem can be solved following arguments similar to those 
presented in section [3] 

The integrated signals (in appropriate units) recorded by experiments 1 and 2 
from a given energy deposition in the atmosphere will be given by 

S ! = E t Y tSl = £iY lSl ; 5 2 = E t Y t e 2 = E 2 Y 2 s 2 . (8) 

In this case Ye includes not only the fluorescence contribution given by equa- 
tion Q but also that of the Cherenkov light. In order to determine the relative 
energy scale of these experiments (E 2 /Ei), the energy of given showers should be 
reconstructed by each experiment using the respective signal. Nevertheless, if we 
assume again that experiment 1 is using the true fluorescence yield, the relation- 
ship 5 2 = £ , iY 1 e 2 holds, that is, we could even calculate the signal of telescope 2 
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from that of telescope 1. Instead of that, the E 2 IEi ratio can be obtained directly 
from expression ([5]) which also applies in this case. While the Ye quantities were 
previously computed analytically for modelled showers (section [3]), now we can 
make the reconstruction algorithm to do it for us using real data, including in this 
way the Cherenkov contribution. 

In general the reconstructed energy is proportional to the reciprocal Ye value, 
the telescope signal being the proportionality constant, i.e., E = S /Ye. From the 
S i signal, E* parameters can be obtained by reconstructing the energy assuming 
in the algorithm the optical efficiency of experiment 2 either with the fluorescence 
yield of experiment 1 or that of experiment 2, 



E\ = —^; E* 2 = . (9) 



_^J_- E* = — 
Yie 2 ' 2 Y 2 e 2 

Note that E\ and E* 2 are not actual values of energy deposition but those of E\ 
and E 2 scaled by the common S1/S2 factor, 

El = ^Eu E* 2 = ^E 2 , (10) 

i3 2 >J 2 

and therefore the relative energy scale can be determined from 

E 2 E* 2 



Ei E\ 



(11) 



From arguments similar to those of section Q, the condition for experiment 
1 of using the true fluorescence yield can be relaxed to just using the true fluo- 
rescence spectrum, although in this case we have to assume that the Cherenkov 
contribution is small with respect to that of fluorescence. 



6 Conclusions 

In this work we have evaluated the effect of the fluorescence yield choice on the 
energy reconstruction of extensive air showers. Using a simple analytical proce- 
dure we found that Auger (without T and h effects) would increase the energy by 
about 9% and 3% if the fluorescence yield of TA and HiRes respectively were 
used instead. 

For the comparison of the energy scales, the optical efficiency of the experi- 
ments also matters. A simple procedure allows to carry out this comparison ana- 
lytically. Assuming that Auger uses the true fluorescence spectrum, i.e., AIRFLY 
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0|, we found that the TA energy is larger than that of Auger (including temper- 
ature and humidity effects) by about 10%. Notice that this energy deviation is 
accidentally lowered because TA does not take into account the temperature and 
humidity contributions to the fluorescence yield. Otherwise the energy scale of 
TA would be 13% larger. 

The above results are somewhat overestimated because the analytical proce- 
dure shown here does not take into account the contribution of Cherenkov light 
in the reconstruction of real data. A procedure has been described that allows to 
make these comparisons from the reconstruction of real data of any experiment. 
Some tests have shown that the analytical AE values are overestimated with re- 
spect to those of real data by about 15% and therefore the energy scale of TA 
would be larger than that of Auger by about 6% (10%) if the humidity and tem- 
perature effects are (not) taken in account in Auger. 

Finally, from the present analysis we also conclude that the energy of TA 
(BRM and LR eyes) is about 10% larger than that of HiRes (i.e., MD eye of 
TA) due to the different fluorescence spectra used by these experiments. 
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